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Abstract 

We consider the dissociation of heavy quarkonium in a medium close to thermal equilibrium 
but with a small momentum space anisotropy. Dissociation is defined to take place when the 
width of the ground state equals its binding energy. We show that if the anisotropic medium 
is obtained isentropically from the equilibrium one, then to first order in the anisotropy pa- 
rameter the dissociation temperature remains unchanged. If, in contrast, the non-equilibrium 
system has a smaller entropy density than the equilibrium one, then the dissociation temper- 
ature increases with respect to the isotropic case, by up to ~ 10% for modest anisotropics. 
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1. The physical picture of quarkonium dissociation within a deconfined medium [T] has 
undergone a slight refinement within the last couple of years. The heavy quark and anti- 
quark are bound together by almost static (off-shell) gluons, and the issue boils down to 
how the gluon self-energy looks at high temperatures. It turns out that the gluon self- 
energy has both a real and an imaginary part. The real part is known as Debye screening, 
while the imaginary part is referred to as Landau damping. Traditionally, it was thought 
that quarkonium dissociates when Debye screening becomes so strong that the corresponding 
Schrodinger equation supports no more bound state solutions. The new suggestion is that 
quarkonium effectively dissociates already at a lower temperature OS], when the binding 
energy is non-zero but overtaken by the Landau-damping induced thermal width [2] . (Finite 
width phenomena have also been discussed in a phenomenological approach [4].) 

More quantitatively, let us denote by M the heavy quark pole mass, by T the temperature, 
by r the quarkonium Bohr-type radius, and by g the QCD coupling constant (a s = g 2 /4ir). 
Then, as is familiar from the hydrogen atom, r ~ l/g 2 M. If quarkonium dissociates through 
Debye screening, the dissociation temperature is parametrically determined by rrriD ~ 1, 
where ttt-d ~ gT is the Debye mass, leading to T ~ gM. However, the imaginary part of the 
real-time static potential [2j overtakes the binding energy already at smaller temperature, 
with mnr < 1; more precisely, in the range g 2 M < T < gM [3]. By formulating the problem 
in an effective theory framework, this estimate was refined to T ~ g 4//3 M in ref. [5], and 
further to T ~ g 4 / 3 ln -1 / 3 (l/g)M in ref. [6j. The last estimate was independently reproduced 
with another approach in ref. [7j. 

To summarise, the Debye-screened picture of quarkonium dissociation seems, in retro- 
spective, overly conservative. To rephrase this more polemically, one can argue that Debye 
screening is not the dominant mechanism responsible for dissociating the quarkonium state 
in a thermal environment [7]. 

All of the considerations above refer to the theoretically transparent situation that the 
medium is in thermal equilibrium. For phenomenological applications it may be interesting 
to study non-equilibrium environments as well. For instance, in ref. [8] the case was consid- 
ered that the quarkonium system has a non-zero velocity with respect to the medium. Here, 
in contrast, we inspect a system in which quarkonium is at rest but the "hard" partonic de- 
grees of freedom have distribution functions in momentum space which contain an anisotropy, 
or a preferred direction. This situation could emerge as a result of Bjorken expansion at the 
early stages of a central heavy ion collision, and it has therefore been extensively studied in 
the literature recently, for instance in the context of plasma instabilities ([9] and references 
therein) as well as for observables such as heavy quark energy loss [10], heavy quark momen- 
tum broadening [TT], photon production [12], dilepton production [13], jet quenching [14], [T5]. 
and, most relevantly for us, the heavy quark potential [16] and solutions of the Schrodinger 
equation [17] . 
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2. Strictly speaking, the suggestion to give the hard modes an anisotropic momentum 
distribution does not lead to a self-consistent non-equilibrium "ground state": rather, the 
system has "tachyonic" modes, meaning that one is trying to compute around a wrong 
extremum. However, as we will see, for our observable these problems are absent in the 
limit of a small anisotropy, and we restrict to this case in the following. 

For a general momentum distribution, the Hard Thermal Loop gluon self-energy obtains 
the form [18] 

d 3 p „„„ „, , /„„ v u k a 



11% [K) = g 2 J ^ v»d a f(p) (V Q 



v ■ K + i0+ ' ( " 
where Minkowskian conventions are assumed, v = (l,p/p), and the subscript R indicates 
that the self-energy appears in the inverse of the retarded propagator. Following ref. |19| . the 
hard mode momentum distribution function is assumed to be of the form 



Hp) = /iso(vy+M) • (2) 

Here £ £ (— 1, oo) is a real parameter, and p z the momentum component in the beam direction; 
the values corresponding to the Bjorken expansion induced anisotropy are £ > 0. 

Carrying out a change of variables to p = v/p 2 + and denoting p z = pcos9 p , the 
momentum integration may be written as 



d 3 
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where v is a unit vector, with v z = cos 6 p ; the angular variables determining the direction of 
v remain unchanged in the substitution; and T is an arbitrary function. 

For the real-time static potential, we need the component of the self-energy, near the 
static limit -C k. Expanding to first order in £ and k° yields 
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where 6^ is the angle between k, z and we denoted 

d 3 p d/iso(p) 



-9 



(2tt) 3 dp 
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3. For the real and imaginary parts of the real-time static potential, we need the gluon 
propagator near the static limit. More precisely, an explicit computation of the static poten- 
tial d la ref. [2] can be rephrased by noting that, if the static limit exists (it does at 0(£) 
but not at 0(£ 2 )), it is the time-ordered gluon propagator that matters [20J, as would be 
expected in the naive real-time formalism: 

_0_k_ e__+e_ 1 iA oo (Q) k) (6) 
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where Cp = {N£ — 1)/2N C and N c = 3. The time-ordered propagator can in turn be written 

as 

iA§? = A°° + 2in s (k°) ImA» « A™ + ^ ImA™ , (7) 

where A^ is the 00-component of the retarded propagator, and ra B is the Bose-Einstein 
distribution function. We have here assumed that, unlike the hard modes, the soft gluons are 
already in thermal equilibrium at a temperature T. It can be verified that at the required 
order in A: and £, the self-energies ng,n^ do not contribute to A^°. Making use of 
eq. ([2]) and expanding to first order in £, we thus obtain that in the static limit 
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Inserting eq. © into eq. ([6]) , and noting that cos Oy. = cos # r cos B^ r + sin 6> r sin cos 0fc r , 
where 9 r is the angle between r, z and 6^, (j)^ are the angular variables between k, r, we can 
integrate over Okn^kr- The potential becomes 
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where r = m^r. Furthermore, we have defined 
p{f) = e 



1 1 1 
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The real part of the potential agrees with the result of ref. [16]. The functions appearing in 
the imaginary part are finite for all r ; <j)\ and <p2 vanish both at small and large distances. 



4. In order to obtain a theoretically consistent result in the weak-coupling regime, it is 
important to systematically account for all the effects to a given order, emerging from the 
various momentum and frequency scales of the problem. This can be achieved by employing 
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effective field theory methods [21] (for a review, see ref. [6]). The upshot is that r is 
formally a small parameter, f < 1 [3], and can be expanded in. 

Expanding in f , the term proportional to £ of the real part of the potential is seen to be 
of 0(g 2 m'^ ) r), i.e. suppressed by two powers of f with respect to the ^-independent term. 
Therefore it can be omitted in the weak-coupling limit. 

In contrast, in the imaginary part (the second line of eq. 0) the corrections from £ are of 
the same order as the leading term: 

f 2 / 1 4 \ 

0(f) « _^i n __ 7E + _j , (14) 

Mr) « y( ln J- 7E + 0' (15) 

f 2 / 1 47\ 

Mf) » y( ln ?-^ + 6o)' (16) 

In particular, at leading logarithmic order, all three functions behave identically. 

5. Before we inspect more precisely the consequences of eq. ([9]), we need to discuss the 
value of the parameter ttid appearing in it, defined by eq. ©. 

Often the function /i so in eq. © is taken to be a sum of Bose-Einstein and Fermi-Dirac 
distribution functions, with a "temperature" T appearing as a parameter (see, e.g., eq. (1) of 
ref. [22J). It is clear, however, that since /i so was part of a specification of a non-equilibrium 
state (cf . eq. ([2]) ) , a precise physical meaning can only be given to T through further arguments 
(cf. ref. [23] and references therein). We have already specified our "late-time" setting above 
(strongly interacting soft modes thermalized at a physical temperature T, weakly interacting 
hard modes still anisotropic), so we need to rethink the meaning of the parameter appearing 
in /i so for this case. We denote this parameter by T' in the following. 

Consider now the entropy density of the system. Making use of the substitution in eq. ([!]) , 
it is given by an integral of the type 

w eq = J Hp) ~ j ^3 Hp) (i - § £ ™s 2 e p ) , (17) 

where we expanded to leading order in £. Carrying out the angular integral, and assuming a 
massless system without chemical potentials, this yields 

Weq = (l - |) CT ' 3 . ( 18 ) 

with a certain constant c. Similarly, the energy density becomes e non -eq = | (1 — 2£/3)cT' 4 . 

We could now envisage two cases. If we define the non-equilibrium state by T 1 = T, then 
we observe that, for £ > 0, it has less entropy density than the corresponding equilibrium 
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state, with s eq = cT 3 . On the other hand, if we impose the physical condition that the 
non-equilibrium state be related "isentropically" to the equilibrium one, then we are left to 
conclude that the parameter T' should be chosen as 



T' 



1 + - )T . 

6/ 



(19) 



The same outcome results if we define the temperature through T = ds non . cq / de non . eq . 



6. In order to estimate the temperature at which quarkonium dissociates, we finally carry 
out a parametric computation according to the discussion above. As we will see the effects of 
the anisotropy parameter £ can be fully accounted for without being concerned about various 
numerical factors. Thus, we estimate the magnitude of the binding energy by the real part 
of the static potential, expanded to leading order in f and evaluated at the distance scale of 
the Bohr radius: 

9 2 C F 



Rc[V> 



4-7rr 



1 + 0{r) 



(20) 



In the imaginary part, in contrast, corrections involving £ are of order unity. Considering the 
s-wave ground state, the term proportional to 1 — 3cos 2 # r in eq. ([9]) averages to zero at first 
order in perturbation theory (corrections will be of order £ 2 ), so to leading- logarithmic order 
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If we now assume the temperature T' parameterizing the non-equilibrium system to be isen- 
tropically obtained from the equilibrium case, leading to eq. (|19p . then the Debye mass 
parameter defined in eq. ([5]) evaluates to 



g 2 T 12 « g z T z ( 1 + 



2rp2 



(22) 



Thus ^-dependence cancels on the right-hand side of eq. (|2ip to C?(£); at leading- logarithmic 
order the dissociation temperature remains at the value 



- melt 
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(23) 



obtained in refs. [HE]. 

If, in contrast, we assume T' 



T, then mj^ ~ 9 T , Temperature appears in eq. ([21 



cubically, meaning that equality with the real part of eq. (|20p is obtained for 
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For £ ~ 1, up to which range our small-£ approach might be assumed qualitatively reasonable, 
we thus obtain a 10% increase in the dissociation temperature. 

7. To conclude, while the precise definition of the heavy quarkonium dissociation tem- 
perature is ambiguous, requiring a convention on the shape of the corresponding smoothly 
evolving spectral function (for results within the weak-coupling expansion, see ref. [Sj); it 
appears that for the class of non-equilibrium states introduced in ref. [H] , the relative change 
caused by a momentum space anisotropy can be estimated analytically at leading-logarithmic 
order. On the other hand, the definition of the non-equilibrium state itself contains a hidden 
ambiguity, in that a temperature-like parameter is introduced which requires further justi- 
fication. We have related this parameter to the entropy density of the system, and thus 
arrived at two physically distinct results, eqs. (|23j) and (|24|) . In general, it is expected that 
non-equilibrium states have less entropy than the equilibrium one, under which conditions 
eq. ([23]) could be a better estimate than eq. ([23]) . indicating an increase of the dissociation 
temperature; nevertheless, reasonable arguments could also be given in favour of eq. (|19p . 
leading to eq. ([23]) . In any case, our results suggest that quarkonium dissociation is primarily 
a probe of the entropy density of the system, i.e. of the number of hard modes, rather than of 
Debye screening, i.e. of soft collective phenomena. At the same time, given the ambiguities 
appearing, it seemed to us that the analytic leading-logarithmic order-of-magnitude estimate 
is about as far as one needs to go for small anisotropies; for larger ones, numerical simulations 
would be needed. 

Note added 

Recently a paper appeared [25] where the same problem is considered as here. The results 
deviate slightly from ours because the physics setting is different: unlike in our eq. ([7]), even 
the soft gluons are assumed to have a non-thermal distribution function, which leads to 
an additional term in the gluon propagator (the 2nd line of eq. (14) in ref. [25]). We note, 
however, that if we re-express within our setting the soft gluon T in terms of T 1 from eq. (I19p . 
and identify T' with the temperature parameter of ref. [25] ; then curiously our eq. (]2T1) does 
reproduce eq. (58) of ref. [25] . 
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